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Galactic Sun's motion in the Cold Dark Matter, Modified Newtonian 
Dynamics and MOdified Gravity scenarios 
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We numerically integrate the equations of motion of the Sun in Galactocentric Cartesian rectangular coordinates for —4.5 
Gyr < t < in Newtonian mechanics with two different models for the Cold Dark Matter (CDM) halo, in Modified 
Newtonian Dynamics (MOND) and in Modified Gravity (MOG) without resorting to CDM. The initial conditions used 
come from the latest kinematical determination of the 3D Sun's motion in the Milky Way (MW) by assuming for the 
rotation speed of the Local Standard of Rest (LSR) the recent value So — 268 km s^^ and the lAU recommended value 
Go = 220 km s"^ the Sun is assumed located at 8.5 kpc from the Galactic Center (GC). For 0o = 268 km s ^ the 
birth of the Sun, 4.5 Gyr ago, would have occurred at large Galactocentric distances (12 — 27 kpc depending on the model 
used), while for 8o = 220 km s^^ it would have occurred at about 8.8 — 9.3 kpc for almost all the models used. The 
integrated trajectories are far from being circular, especially for Oo = 268 km s~^, and differ each other with the CDM 
models yielding the widest spatial extensions for the Sun's orbital path. 
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1 Introduction 

In several astrophysical systems like, e.g., spiral galaxies 
and clusters of galaxies a discrepancy between the observed 
kinematics of their exterior parts and the predicted one on 
the basis of the Newtonian dynamics and the matter de- 
tected from the emitted electromagnetic radiation (visible 
stars and gas clouds) was present since the pioneering stud- 
ies by Zwicky (1933) (he postulated the existence of un- 
detected, baryonic matter; today, it is believed that the hid- 
den mass is constituted by non-baryonic, weakly interacting 
particles) on the Coma cluster, and by Bosma (1981) and 
Rubin et al. (1982) on spiral galaxies. More precisely, such 
an effect shows up in the galactic velocity rotation curves 
(Persic & Salucci 1996a, 1996b) whose typical pattern af- 
ter a few kpc from the center differs from the Keplerian 
1/y/r fall-off expected from the usual dynamics applied to 
the electromagnetically-observed matter. 

As a possible solution of this puzzle, the existence of 
non-baryonic, weakly-interacting Cold Dark (in the sense 
that its existence is indirectly inferred only from its gravita- 
tional action, not from emitted electromagnetic radiation) 
Matter (CDM) was proposed to reconcile the predictions 
with the observations (Rubin 1983) in the framework of the 
standard gravitational physics. 

Oppositely, it was postulated that the Newtonian laws 
of gravitation have to be modified on certain acceleration 
scales to correctly account for the observed anomalous kine- 
matics of such astrophysical systems without resorting to 
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Still undetected exotic forms of matter. One of the most 
phenomenologically successful modification of the inverse- 
square Newtonian law, mainly with respect to spiral galax- 
ies, is the Modified Newtonian Dynamics (MOND) (Mil- 
grom 1983a, 1983b, 1983c) which, as we will see below, 
yields aw l/r gravitational force for very small acceler- 
ations. 

Another alternative theoretical framework recently put 
forth to also explain, among other astrophysical and cos- 
mological observations, the observed kinematics of the out- 
skirts of galaxies without resorting to CDM is MOdified 
Gravity (MOG) (Moffat & Toth 2008) which yields a long- 
range Yukawa-like modification of the Newtonian inverse- 
square law. 

In this paper we want to investigate the motion of the 
Sun in the Milky Way (MW) over timescales of the order of 
its lifetime, i.e. 4.5 Gyr, in CDM, MOND and MOG scenar- 
ios in view of the recent accurate measurements of the solar 
kinematical parameters (Reid et al. 2009). This may help in 
discriminating, at least in principle, the different theoreti- 
cal scenarios examined; moreover, it will be interesting to 
see if the Galactic orbital motions of the Sun backward in 
time predicted by CDM, MOND and MOG are compatible 
with known and accepted knowledge concerning its forma- 
tion and the evolution of the life in our planet. This may 
also yield an independent consistency test of the new kine- 
matical determinations of MW (Reid et al. 2009). Finally, 
we note that the approach presented here may be used, in 
principle, also for other dynamical features. 
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Table 1 Parameters of the disk and bulge models by Xue 
et al. (2008). 



Table 3 Parameters of the logarithmic CDM halo model; 
values by Law et al. (2005). 
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Table 2 Parameters of the CDM NFW halo model by Xue 
et al. (2008). The values quoted for r^ii and c come from an 
average of those by Xue et al. (2008). 
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2 The models used 
2.1 The CDM NFW model 

The CDM model tested by Xue et al. (2008) with several 
Blue Horizontal-Branch (BHB) halo stars consists of three 
components. Two of them are for the disk 

GA/disk [l-exp(-i)] 

Udisk = , (1) 

r 

where b is the disk scale length, and the bulge 

J J _ GA/buigc 



r + co 



(2) 



where c is the bulge scale radius. The third component is for 
the CDM NFW halo (NavaiTO et al. 1996) 
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NFW 



AtKJpsT^^^ / cr 



with Tvir is the radius parameter and 



3 ln(l + c)-^' 



(3) 



(4) 



in which Jim is the fraction of matter (including baryons and 
DM) to the critical density, 6th is critical overdensity of the 
virialized system, c is the concentration parameter, and 



Per = 



3Hl 
8nG 



(5) 



is the critical density of the Universe determined by the 
Hubble parameter Hq at redshift 0. The values used for the 
parameters entering eq. ([T]i and eq. (|2|i are in Table [T] those 
entering eq. (O are in Table|2] 

2.2 The logarithmic CDM halo 

Another model widely used consists of the Miyamoto & Na- 
gai (1975) disk 



disk 



(,GM, 



disk 



(6) 



where k and B are the disk scale length and the disk scale 
height, respectively, the Plummer (1911) bulge 



^bulgc 



GMbuigt 
r + C 



(7) 



where C is the bulge scale radius, and the logarithmic CDM 
halo (Binney & Tremaine 1987) 



(f) , (spherical halo) (8) 



where d is the DM halo scale length and Whaio describes 
the DM halo dispersion speed (which is related to the total 
DM halo mass). The parameters's values are in Table[3] The 
masses of the disk and the bulge used by Law et al. (2005) 
are those by Johnston et al. (1999): A/disk = 1 x 10" M©, 
A^buigo — 3.4 X lO^*' M0 yielding a total baryonic mass 
of Af = 1.34 X 10^^ M0; however, such a value is almost 
twice the most recent estimate (Af = 6.5 x 10^" Mq) by 
McGaugh (2008) who includes the gas mass as well and 
yield Afdisk = 2.89 x 10^" Mq and Afbuigo = 2.07 x 10^" 
Mq. The model of eq. ©-eq. ([8]), with the parameters' val- 
ues by Law et al. (2005) and Johnston et al. (1999), has been 
recently used by Willett et al. (2009) to study the motion 
of the Grillmair & Dionatos (2006) tidal stellar stream at 
galactocentric distance of r « 16 — 18 kpc; Read & Moore 
(2005) used it to study the motion of the tidal debris of the 
Sagittarius dwarf at 17.4 kpc from the center of MW. More 
specifically, the CDM halo model of eq. ^ corresponds to 
a CDM halo mass 



halo 



2^Lio^ 



G(r2 +d2)' 



(9) 



so that 

A//halo(r - 60 kpc) ^ 3.5 X 10" Mq, (10) 

in agreement with the value by Xue et al. (2008) 

A4aio(?- = 60 kpc) = (4.0 ± 0.7) x lO" Mq. (11) 

2.3 Modified Newtonian Dynamics 

MOND postulates that for systems experiencing total grav- 
itational acceleration A < Aq, with (Begeman et al. 1991) 



Aq = (1.2 ± 0.27) X 10"^° m 



A Amond 
More precisely, it holds 
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fi{X) — > 1 for X ^ 1, i.e. for large accelerations (with 
respect to Aq), while fJ.{X) X yielding eq. ( fT3]) for X ^ 
1, i.e. for small accelerations. Concerning the interpolating 
function i.J.{X), it recently turned out that the simple form 
(Famaey & Binney 2005) 



X 



1 + X 



(15) 



yields very good results in fitting the terminal velocity curve 
of MW, the rotation curve of the standard external galaxy 
NGC 3198 (Famaey & Binney 2005, Zhao & Famaey 2006, 
Famaey et al. 2007) and of a sample of 17 high surface 
brightness, early-type disc galaxies (Sanders & Noorder- 
meer 2007); eq. ( fT4l l becomes 



^Newton M 



(16) 



2 \ V ^Newton / 

with eq. ( fTSl l. Thus, in the following we will use eq. (fTST l. 
2.4 Modified Gravity 

MOG is a fully covariant theory of gravity which is based on 
the existence of a massive vector field coupled universally 
to matter. The theory yields a Yukawa-like modification of 
gravity with three constants which, in the most general case, 
are running; they are present in the theory's action as scalar 
fields which represent the gravitational constant, the vector 
field coupling constant, and the vector field mass. An ap- 
proximate solution of the MOG field equations by Moffat 
& Toth (2009) allows to compute their values as functions 
of the source's mass. The resulting Yukawa-type modifica- 
tion of the inverse-square Newton's law in the gravitational 
field of a central mass M is 



A = - 



{l + a[l-(l + Hexp(-H]}, (17) 



with 
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c' 
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(18) 



(19) 



where Gn is the Newtonian gravitational constant and 

Goo ^ 20Gn, (20) 
C[ ~ 25,OOO0Z^, (21) 



G, 



6,250 




(22) 



Such values have been obtained by Moffat & Toth (2009) 
as a result of the fit of the velocity rotation curves of some 
galaxies in the framework of the searches for an explana- 
tion of the rotation curves of galaxies without resorting to 
CDM. The validity of eq. ( fTTI l in the Solar System has been 
recently questioned by lorio (2008). 

For M = 6.5 X 10^° Mq we have the values of TableH 



Table 4 MOG parameters a and A for M = 6.5 x 10^ 

Mq. 
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(kpc) 



16 



41 



Table 5 Galactocentric initial conditions for the Sun 
(Reid et al. 2009); the positive y axis is directed from the 
Galactic Center (GC) to the Sun, the positive x axis is di- 
rected toward the Galactic rotation, the positive z axis is 
directed toward the North Galactic Pole (NGP). We have 
used (Reid et al. 2009) 8o = 268 km s^^ for the rotation 
speed of LSR, and Uq = 10.3 km s^^, Vq = 15.3 km s-\ 
Wq — 7.7 km s^^ for the standard solar motions toward 
GC, ^ = 90 deg and NGP, respectively; thus, with our con- 
ventions, ±0 = ^0 + Qo, 2/0 = —Uo, zq = Wq. See Fig. 7 
of (Reid et al. 2009). 
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3 Solar motions backward in time 

We will numerically integrate witlQ MATHEMATICA the 
Sun's equations of motion in Cartesian rectangular coordi- 
nates for —4.5 Gyr < t < wit the initial conditions of 
Table |5l in it we use the recently estimated value Qq = 268 
km s^^ for the rotation speed of the Local Standard of Rest 
(LSR) (Reid et al. 2009). 

We will use both the NFW and logarithmic halo models 
for CDM which have been tested independently of the solar 
motion itself. In particular, we will numerically integrate the 
three scalar differential equations, written in cartesian coor- 
dinates, corresponding to the vector differential equation 



(23) 



where U is the sum of eq. ([T]i-eq. ^ with the values of Ta- 
ble |2] in the first case (NFW), and of eq. ©-eq. dHJ with 
the values of Table |3] in the second case (logarithmic); we 
use the figures of Table [T] in both cases. For the two mod- 
els of modified gravity considered we will put the carte- 
sian components of eq. (MOND) and eq. {TT^ (MOG) 
in the right-hand-sides of the differential equations to be 
integrated. Concerning MOND, let us note that Eq. (fl4] l. 
from which eq. (fTSl l has been derived by means of eq. (fTsT i 
for fJ.{X), strictly holds for co-planar, spherically and axi- 
ally symmetric mass distributions; otherwise, the full mod- 
ified (non-relativistic) Poisson equation (Bekenstein & Mil- 
grom 1984) 



= 47rGp 



(24) 



We will use the default options of NDSolve. 
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y, kpc 




Fig. 1 Section in the {x, y} plane of the numerically inte- 
grated trajectories of the Sun for a) CDM logarithmic halo 
(dash-dotted blue line) b) CDM NFW halo (dotted yellow 
hne) c) MOND with ^ = X/{l + X) (dashed green line) c) 
MOG (dashed red line). The time span is —4.5 Gyr < t < 0. 
The initial conditions are those of Table |5] with (Reid et 
al. 2009) Go = 268 km s^^. 



z, kpc 




Fig. 2 Section in the {x, z} plane of the numerically inte- 
grated trajectories of the Sun for a) CDM logarithmic halo 
(dash-dotted blue line) b) CDM NFW halo (dotted yellow 
line) c) MOND with X/(l + X) (dashed green line) c) 
MOG (dashed red line). The time span is -4.5 Gyr <t<0. 
The initial conditions are those of Table |5] with (Reid et 
al. 2009) eo = 268kms-i. 



must be used. However, in this case, the baryonic mass dis- 
tributions of the bulge and the disk (see eq. ([TJ-eq. (|2]i and 
eq. ©-eq. (I7|) just satisfy the symmetry conditions that al- 
low to use eq. (fT4l) . Moreover, An Mo ~ 1 for M 6, 5 x 
10^" M0. This allows to neglect the so-called External Field 
Effect (FEE) (Sanders & McGaugh 2002, Milgrom 2008) 
because it should amount to about A^/Ao = 0.01 — 0.02 
for MW (Wu et al. 2008). Indeed, it can be shown that, for 
eq. ( fTSt . the MOND acceleration, including FEE, is 



An 
2 



An 



A. 

An 



4Ao 
An 



1 + 



(25) 



which just reduces to eq. ( fTSI ) when A^ ^ ^n, Ac <C Aq, 
as in this case. 

In Figure [T]- Figure [5] we plot the orbital sections in the 
Galactocentric coordinate planes of the numerically inte- 
grated trajectories of the Sun from now to 4.5 Gyr ago. The 
orbital trajectories of the Sun are quite different in CDM, 
MOND and MOG; the largest spatial extension occur for 
the CDM models, especially for the one characterized by 
the CDM logarithmic halo The Galactocentric coordinates 
and distances of the Sun at t = —4.5 Gyr are quoted in 
Table |6l It can be noted that the models considered place 
the birth of the Sun from a minimum of 12.0 kpc (MOND) 
to a maximum of 25 — 27 kpc (CDM) from GC. This may 
pose problems concerning the birth of the Earth itself and 
the development of complex life on it because of, e.g., the 
presumable low level of metallicity (the amount of elements 
heavier than hydrogen and helium) in so distant regions. For 
a discussion and the implications of the concept of Galactic 
Habitable Zone (GHZ), see Lineweaver et al. (2004), Go- 
nalez (2005) and Prantzos (2008). In particular, see Figure 
5 by Prantzos (2008) which shows that the probability of 
having Earths formed 4 and 8 Gyr after the formation of 
MW, which roughly corresponds to the birth of the Sun by 
assuming a Milky Way age of about 10 — 12 Gyr, is prac- 




y, kpc 



Fig. 3 Section in the {y, z} plane of the numerically inte- 
grated trajectories of the Sun for a) CDM logarithmic halo 
(dash-dotted blue line) b) CDM NFW halo (dotted yellow 
line) c) MOND with /i = X/(l + X) (dashed green line) c) 
MOG (dashed red line). The time span is —4.5 Gyr < t < 0. 
The initial conditions are those of Table |5] with (Reid et 
al. 2009) eo = 268kms-i. 



Table 6 Galactocentric coordinates and distance of the 
Sun at < = —4.5 Gyr according to our numerical integra- 
tions of the equations of motion in CDM (logarithmic halo), 
CDM (NFW halo), MOND and MOG for (Reid et al. 2009) 
8o = 268 km s^^. The initial conditions of Table |5] have 
been used. 
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y, kpc 




Fig. 4 Section in the {x, y} plane of the numerically inte- 
grated trajectories of the Sun for a) CDM logarithmic halo 
(dash-dotted blue line) b) CDM NFW halo (dotted yellow 
hne) c) MOND with ^ = X/{l + X) (dashed green line) c) 
MOG (dashed red line). The time span is —4.5 Gyr < t < 0. 
The initial conditions are those of Table |5] with the lAU rec- 
ommended value 8o ~ 220 km s^^. 



z, kpc 




Fig. 5 Section in the {x, z} plane of the numerically inte- 
grated trajectories of the Sun for a) CDM logarithmic halo 
(dash-dotted blue line) b) CDM NFW halo (dotted yellow 
Hne) c) MOND with X/{l + X) (dashed green line) c) 
MOG (dashed red line). The time span is —4.5 Gyr < t < 0. 
The initial conditions are those of Table |5] with the lAU rec- 
ommended value 8o = 220 km s^^. 

tically null at Galactocentric distances larger than 10—15 
kpc. 

In Figure |4l-Figure|6] we repeat the same integrations by 
using the lAU recommended value Oq = 220 km s^^. The 
Galactocentric coordinates and distances of the Sun at < = 
-4.5 Gyr for Oo = 220 km s^^ ai-e quoted in Table |3 
Now the situation is quite different because, apart from the 
logarithmic CDM halo, all the other models locate the birth 
of the Sun at 8.8 - 9.3 kpc from GC 

4 Summary and Conclusions 

By using the latest kinematical determinations of the full 
3D solar motion in the Milky Way, implying a LSR rotation 



z, kpc 




Fig. 6 Section in the {y, z} plane of the numerically inte- 
grated trajectories of the Sun for a) CDM logarithmic halo 
(dash-dotted blue line) b) CDM NFW halo (dotted yellow 
line) c) MOND with ^ = X/{l + X) (dashed green line) c) 
MOG (dashed red line). The time span is -4.5 Gyr <t<0. 
The initial conditions are those of Table|5]with the lAU rec- 
ommended value 6o — 220 km s^^. 

Table 7 Galactocentric coordinates and distance of the 
Sun at t = —4.5 Gyr according to our numerical integra- 
tions of the equations of motion in CDM (logarithmic halo), 
CDM (NFW halo), MOND and MOG. The initial condi- 
tions of Table |5]have been used with the lAU recommended 
value Go = 220 km s'^. 
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-18.1 


-4.7 


0.1 


18.8 


CDM NFW 
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0.2 
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MOND (m=t^) 
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0.2 


8.9 
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9.2 


1.1 
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9.3 



of &0 = 268 km s~^, we used them as initial conditions 
for numerically integrating the equations of motion of the 
Sun backward in time (—4.5 Gyr < t < 0) in Newtonian 
mechanics with two models for the CDM halo, in MOND 
and in MOG. As a result, the orbital trajectories are not cir- 
cular and differ each other, with the CDM models yielding 
the widest spatial extensions of the Sun's orbit. It turns out 
that for t = —4.5 Gyr the Sun is at quite large Galactocen- 
tric distances (r « 12 — 27 kpc depending on the models). 
Instead, by using the standard lAU value 8o = 220 km 
s^^ the situation is different: the orbits are less wide and at 
t = -4.5 Gyr the Sun is at 8.8- 9.3 kpc from GC for almost 
all the models considered. 
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